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Abstract 

We find new classes of exact solutions to the Einstein-Maxwell system of equa- 
tions for a charged sphere with a particular choice of the electric field intensity 
and one of the gravitational potentials. The condition of pressure isotropy is 
^Jy reduced to a linear, second order differential equation which can be solved in 

general. Consequently we can find exact solutions to the Einstein-Maxwell field 
equations corresponding to a static spherically symmetric gravitational potential 
in terms of hyper geometric functions. It is possible to find exact solutions which 
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can be written explicitly in terms of elementary functions, namely polynomials 
and product of polynomials and algebraic functions. Uncharged solutions are 
regainable with our choice of electric field intensity; in particular we generate 
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the Einstein universe for particular parameter values. 
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1 Introduction 

In recent years a number of authors have found solutions to the Einstein-Maxwell field 
equations for static spherically symmetric gravitational fields with isotropic matter. 
These exact solutions must match at the boundary to the unique Reissner-Nordstrom 
metric which is the exterior spacetime for a spherically symmetric charged distribution 
of matter. The models generated are used to describe relativistic spheres with strong 
gravitational fields as is the case in neutron stars. It is for this reason that many inves- 
tigators use a variety of techniques to attain exact solutions. A comprehensive list of 
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Einstein-Maxwell solutions, satisfying a variety of criteria for physical admissability, is 



provided by Ivanov [1]. The exact solutions may be used to study the physical features 
of charged spheroidal stars as demonstrated by Komathiraj and Maharaj |2y, Sharma 



et al |3], Patel and Koppar [4J, Patel et al [5j, Tikekar and Singh [6j and Gupta and 
Kumar [7|. These analyses indicate that the Einstein-Maxwell exact solutions found 
are relevant to the description of dense astronomical objects. Some other individual 
treatments include the Sharma et al [s| study of cold compact objects, the Sharma 
and Mukherjee 0] consideration of strange matter, and the Sharma and Mukherjee 



lOj analysis of quark-diquark mixtures in equilibrium. Charged relativistic spheres 



may be used to model core-envelope stellar configuration as shown by Thomas et al 



11], Tikekar and Thomas [12|, and Paul and Tikekar [13J where the core consists an 



isotropic fluid and the envelope comprises an anisotropic fluid. 

In order to integrate the field equations, various restrictions have been placed by 
investigators on the geometry of spacetime and the matter content. Mainly two distinct 
procedures have been adopted to solve these equations for spherically symmetric and 
static manifolds. Firstly, the coupled differential equations are solved by computation 
after choosing an equation of state. Secondly, the exact Einstein-Maxwell solution can 
be obtained by specifying the geometry and the form of the electromagnetic field. We 
follow the latter technique in an attempt to find solutions in term of special functions 
and elementary functions that are suitable for the description of relativistic charged 
stars. This approach was first used by John and Maharaj |l4j that yielded an uncharged 
star which approximates a polytrope close to the centre. This particular exact solution 
was extended to a wider class of solutions by Maharaj and Thirukkanesh 15| in the 
presence of charge. Also Thirukkanesh and Maharaj 16| found a family of Einstein- 
Maxwell solutions that contain the Durgapal and Bannerji 17J neutron star model. 



Komathiraj and Maharaj [2| presented a general class of Einstein-Maxwell solutions 
that contain Tikekar [3] spheroidal stars as a special case which are physically viable 
neutron star models. Hence the approach followed in this paper has proved to be a 
fruitful avenue for generating new exact solutions for describing the interior spacetimes 
of charged spheres. 

The objective of this paper is to provide systematically a rich family of Einstein- 

I I 

Maxwell solutions similar to the recent treatment of Komathiraj and Maharaj |2j. 
In Section 2, we rewrite the Einstein-Maxwell equations as a new set of differential 
equations utilising a transformation due to Durgapal and Bannerji [171 ] . We choose 
particular forms for one of the gravitational potentials and the electric field inten- 
sity, which enables us to obtain the condition of pressure isotropy in the remaining 
gravitational potential in Section 3. This is the master equation which determines 
the integr ability of the system. In Section 4, we integrate the condition of pressure 
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isotropy, for particular parameter values, and consequently produce Einstein-Maxwell 
solutions in terms of elementary functions. We demonstrate that exact solutions to the 
Einstein-Maxwell system in terms of hypergeometric functions are possible in Section 
5. In Section 6, we generate two linearly independent classes of solutions by deter- 
mining the specific restriction on the parameters for a terminating series; the general 
solution can be written explicitly in terms of elementary functions. We demonstrate 
that uncharged solutions are regained in the appropriate limit. In Section 7 we discuss 
the physical features, and plot the gravitational and matter variables to show that the 
model is physically acceptable. Finally in Section 8, we show that other solutions, 
outside the class considered in this paper, exist to the Einstein-Maxwell system. 



2 Field equations 

We assume that the interior of a dense compact relativistic star should be spherically 
symmetric. Therefore there exists coordinates (t, r, 6, <p) such that the line element is 
of the form 

fa* = _ e 2u(r) dt 2 + e 2X(r) dr 2 + ^2 + ^2 



The Einstein-Maxwell field equations govern the behaviour of the gravitational field in 
the presence of an electromagnetic field. The Einstein-Maxwell system becomes 



1 , 1( 2A' 

r 2 v r 2 



-(l-e- 2A ) + — e" 2A = p + -E 2 (2a) 



4(1-^ 2A ) + — e- 2X = V -\e 2 ( 2b ) 
r r 2 

e -2A(^ + ^ + 3-z/A'-^) = p+ l -E 2 (2c) 

1 



r 



^e-\r 2 E)' (2d) 



for the line element ([T]). The energy density p and the pressure p are measured relative 
to the comoving fluid 4- velocity u a = c~ u 5q, E is the electric field intensity, a is the 
proper charge density, and primes denote differentiation with respect to r. We are 
utilising units where the coupling constant = 1 and the speed of light c = 1. 

A different but equivalent form of the field equations is generated if we introduce 
new variables 

A 2 y 2 (x) = e 2u{r \ Z(x) = e - 2A(r) , x = Cr 2 (3) 
where A and C are arbitrary constants. Under the transformation ([3]) due to Durgapal 
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and Bannerji [17j, the system (j2J) becomes 

1^-22 = £ + ^ (4a) 
x C 2C v ; 

« Z - 1 p -E 2 
4Z y + — = §"2C (4t) 

4Zx 2 y + 2Za; 2 |/ + ( Za; - Z + 1 - — )y = (4c) 

£ = f (*i? + £) 2 (4d) 

where dots denote differentiation with respect to x. The system of equations (j4j) 
determines the behaviour of gravity for a charged perfect fluid. When E = we 
regain Einsteins equations for a neutral fluid. In the above we have a system of four 
equations in the six unknowns p,p, E, a, y and Z. We are free to specify two of the six 
unknowns; in this treatment we assume forms for Z and E. Once the metric function 
Z and the electric field intensity E are specified then the metric function y can be 
found by integrating Pel) which is then second order and linear in y. The remaining 
unknowns are then obtained from the rest of the system. This is the approach that we 
follow in this paper. Hence the differential equation pel) is the master equation whose 
integration is necessary to determine an exact solution. 



3 Master equation 

We study a particular form of the Einstein-Maxwell system (T4j) by making explicit 
choices for Z and E. For the metric function Z we make the choice 

where k is a real constant. Note that the choice (jSJ) ensures that the metric function 
e 2A is regular and finite at the centre of the sphere. When k = 1, in the absence of 
charge, we regain the Schwarzschild interior metric. Also observe that when k = we 
regain the metric function considered by Hansraj and Maharaj [3] which generalises 
the Finch and Skea j2o| neutron star model. We have chosen the form (jSJ) as it provides 

r L 

for a wider range of possibilities than the solutions of Hansraj and Maharaj [19J, and 
it does produces charged and uncharged solutions which are necessary for a realistic 
model. 

On substituting dSJ) in pel) we obtain 



4(1 + A;a;) 2 (l + x)y + 2(1 + kx)(2k - 1 + kx)y + 



[1-kf- 



E 2 (l + x) 2 
~C~x 



y = (6) 
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It is convenient at this point to introduce the following transformation 

± + x = KX, ^ = K, y(x) = Y(X) (7) 

This transformation enables us to rewrite the second order differential equation ([6]) in 
a simpler form. Under the transformation (JTj), equation ([6]) becomes 



4X%X _ ,)iL + 2X(X _ 2) g 



K 



E 2 K(K + 1) 2 (X 
C[K(X - 1) - 



y = o 



in terms of the new dependent and independent variables Y and X respectively. 

It is necessary to specify the electric field intensity E to integrate (jSJ). A variety of 
choices for E is possible but only a few are physically reasonable which generate closed 
form solutions. We can reduce (|SJ) to simpler form if we let 

E 2 a[K(X - 1) - 1] _ olKx 

~C ~ K(K + 1) 2 (X - l) 2 ~ (K + 1) 2 (1 + x) 2 ^> 

where a is a constant. The form E 2 in (Q is physically palatable because E remains 
regular and continuous throughout the sphere. In addition the field intensity E vanishes 
at the stellar centre, and has positive values in the interior of the star for relevant choices 
of the constants a and K. Upon substituting the choice ([9]) in equation (jSJ) we obtain 

d 2 Y dY 

4x2{x ~ l) lx 2 + 2X{X ~ 2) dx + {K ~ a)Y = (10) 

which is the master equation for the system When a = there is no charge. 
Equation ([TO]) has to be integrated to find an exact model for a charged sphere. 



4 Special case : elementary functions 



We can immediately integrate (TIOl) for the special case K = a ^ 0. Equation ( jTOl) is 
separable and we obtain the solution 



Y(X) = Cl {VX^T - arctan VX~^1) 



c 2 



where C\ and c 2 are constants of integration. In terms of the independent variable x 
we can write 

1 ~\~ x 1 1 ~t~ x 

y(x) =cx[\/ — arctan J ) + c 2 
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Hence the complete solution of the Einstein-Maxwell system (j3J) is then given by 
(K+l) 2 (l + x) 



-,2A 



'2v 



P_ 

c 

v_ 

c 



£2 
~c 



+ 



(K + l + x) 



A 2 



Ci 



1 + x 1 1 + X 

arctan 



K 



K 



+ c 2 



K 2 (6 + x) -6(1 + a:) 



2(K + 1) 2 (1 + x) 2 



2d(K + 1 + x) 



VK(K + l) 2 yTT^ ci */i±* - arctan */*±* + 



c 2 



2(1 + x) 2 -K 2 (2 + x) 
2(K+ 1) 2 (1 + x) 2 

(X+ 1)2(1 +x) 2 



11a) 



lib) 



:iici 



(nd) 

(lie) 



Note that the charged solution (TIT!) does not have an uncharged analogue as the electric 
field intensity E cannot vanish (expect at the centre). This effect essentially results 
from our condition that a = 0). This means that this solution models a sphere 
that is always charged and hence cannot attain a neutral state. A particular class in 
the family of solutions found by Hansraj and Maharaj [19] also demonstrates the same 
feature and E ^ 0. The model (TTTj) is a simple solution of the Einstein-Maxwell system 
which is expressed in terms of elementary functions. 

5 General case : series solution 

With a 7^ K, equation (ITU]) is difficult to solve. However it can be transformed to 
a hypergeometric differential equation which can be integrated using the method of 
Frobenius. We now introduce a new function U(X) such that 

Y(X) = X a U(X) (12) 

where a is a constant. On substituting ffl2l in ffTOl) we obtain 

d 2 U dU 
4X 2 (X-1)—— +2X[(4a+l)X-2(2a+l)] — + [2a(2a-l)X+K-a-Aa 2 ]U = (13) 
d X dX 

We observe that there is considerable simplification if we make the choice 

K - a = 4a 2 (14) 
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This then gives 

d 2 U dU 
2X{X - I)— + [(4a + l)X - 2(2a + 1)]— + a(2a - l)U = (15) 
dX dX 

which is a second order differential equation in terms of the new dependent variable 
U and independent variable X. When a = then a = K and we regain the result of 
Section 4. Therefore we take a ^ in this section to ensure that a ^ K. 
If we let z = 1 — X then (1151) becomes 



z(l-z) 



d 2 U 
~dz* 



1\ 1 

2a + - 2 + - 
2/ 2 



dC/ / 1 . ,. 



(16) 



The result ( TT6l) is a special case of the hypergeometric equation which can be solved 
explicitly in terms of special functions JJ\ and t/ 2 . These special functions are hyper- 
geometric functions and are given by 



and 



T, = /•' [a, a - -, i, : 



U 2 = z^F (a + l,a + l b -,z) 



(17) 



It is now possible to write the solution of ( fl6l) explicitly as a series using the definitions 
of (TlTl) and ( fl8l) . With the help of (JTj) and (fl2l) we obtain the expressions 



X 



if + l + x\ a 

oo 

i=l p=l 



[(p - l)(2p + 4a - 3) + a(2a - 1)] fl + x 
p{2p-3) 



(19) 



and 



V2 [X) 



K+l + x\ a fl + 



x 



K J \ K 

[(2p + l)(p + 2a) + a(2a - 1)] fl + x 

i=i 



i+E(-D-ii 

p=i 



p(2p + 3) 



(20) 



as linearly independent solutions of ([6]). Thus the general solution to the differential 
equation (jSJ), for the choice of the electric field is given by 



y(x) = Aiyi(ac) + A 2 2/ 2 (;r) 

where Ai and A 2 are arbitrary constants, K = -^p, 
by ffT9l) and fl20|) respectively 



if -a 
4 ' 



(21) 

and are y 2 are given 
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From (I2ip and (j4j) we can write the exact solution of the Einstein-Maxwell system 
in the form 

e - (K + i + xy (22a) 

e 2v = A 2 y 2 (22b) 
p [K 2 - 1) (3 + x) - x(5 + 3x) aKx 



C (K+ 1) 2 (1 + x) 2 2{K + 1) 2 (1 + x) 2 

p A(K + 1 + x) 2 y 1 - K 2 + x aKx 



(22c) 



■,:■>,-, ' .., ; .w.- ■,;■>;■, ( 22d ) 

(22e) 



C (K + 1)2(1 + x)y {K + 1) 2 (1 + x) 2(K + 1) 2 (1 + x) 
E 2 aKx 



C (K + 1) 2 (1 + xf 

Unlike the solution presented in Section 4, the models found in this section cannot be 
written in terms of elementary functions in general as the series in (1171) and ( TT81) do 
not terminate. However terminating series are possible for particular values of a, which 
leads to elementary functions, as we show in Section 6. 



6 Elementary functions 

The general solution (|2"T|) is given in the form of a series and can be expressed in 
terms of hypergeometric functions which are special functions. It is well known that 
hypergeometric functions can be written in terms of elementary functions for particular 
parameter values. This statement is also true for the solution found in Section 5 for 
particular values of the parameter a as the two series terminate. Consequently two 
sets of general solutions in terms of elementary functions can be found by restricting 
the range of values of a so that the series terminates. The elementary functions, found 
in this way, are expressible as polynomials and product of polynomials with algebraic 
functions. We can express the first category of solutions, in terms of the original 
variable 

»(,) = Aj K V* (-!)"(» -1) fl+*- 



K+l + xJ ^ (2i)!(2n-2i + l)! \ K 



n-l 



2 \K+1+x)\k) ^ (2i + 3)!(2n-2i-2)! V K ) K ' 
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where K 



a 



An 2 . The second category of solutions is given by 



y{x) 



K 



+ A, 



K + 1 + x 
K 

K+l + x 
4n(n - 1) + 1 



E 



iy- 1 (2i- 1) fl+x 



i=0 

1 + x 



(2i)!(2n - 2% 

| n-2 

E 



K 



K 

-l)\i + 1) 



i=0 



1 + x 



(2i + 3)!(2n-2i-3)! \ # 



(24) 



where K — a 

Therefore two categories of solutions in terms of elementary functions can be ex- 
tracted from the general series in Section 5. The solutions in ( 123]) and (124"]) have a 
simple form, and they have been expressed completely as combinations of polynomials 
and algebraic functions. This has the advantage of simplifying the investigation into 
the physical properties of a dense charged star. As the metric function (j3J) and the 
electric field intensity (Q have not been considered before, we believe that the Einstein- 
Maxwell solutions found here have not been published previously. It is interesting to 
observe that our treatment has brought together the charged and uncharged models 
for a relativistic star. If we set a = in the Einstein-Maxwell solutions (l23j) and (1241) 
then we obtain the solutions for the uncharged case directly. Thus our approach has 
the welcome feature of producing uncharged solutions when E — 0; it is possible that 
the uncharged solutions produced in this procedure may be new. 

We illustrate this feature with an example. We observe that when K—a = 4(n = 1), 
(1231) becomes 

ai {K + 3 + 3x) +oa(l +x)§ 

y{x) = kTTTx (25) 

where a% and 02 are constants. On substituting (|25|) in (|22|) we obtain the general 
solution to the Einstein-Maxwell system of equations as 

(K + l) 2 (l+x) 



,,2A 



■2v 



c 

p_ 

c 



A 



(K+l + x) 2 

ai(K + 3 + 3x) + a 2 (l + x) 



K + l+x 



+ 



£2 
~C 



6(K 2 - 1) + x[(K + 6)(K - 2) - 6a] 
2{K+ 1) 2 (1 + x) 2 

2(K +l + x) [A ai K + a 2 VT z \ r x{3K + 1 + x) 
(K + 1) 2 (1 + x)[ ai {K + 3 + 3z) + a 2 (l + x)l 

2(1 + x)(l - K 2 + x) + K(K - A)x 
2{K+ 1)2(1 + x ) 2 

K(K - A)x 



(26a) 



(26b) 



(26c) 



(K + 1) 2 (1 + x) 



(26d) 
(26e) 
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for our chosen parameter values. When a = 0(K = 4) the electromagnetic field 
vanishes and we get 



,2A 
P_ 

c 

p_ 

c 



25(1 + x) 
(5 + x) 2 

ai(7 + 3x) + a 2 (l + x) 



A 



5 + x 

45 + x(10 - 3x) 
25(1 + x) 2 



ai[3z(x - 2) + 55] + a 2V ^Tx[x(22 + 3x) + 115] 
25(l + x)[ai(7 + 3a;) + a 2 (l + x)%] 



(27a) 
(27b) 

(27c) 
(27d) 



Thus we have generated the uncharged solution (j2"Tj) from the charged solution (T26l) . 



7 Physical features 

We make some brief comments relating to the physics found in this paper. In the general 
solution ([221 . when studying models of charged spheres, we should consider only those 
values of K for which the the energy density p, the pressure p and the electric field 
intensity E are positive. Our choice of the gravitational potential (j5J) is clearly positive 
for a wide range of the parameter values of K. Since y(x) = Aiyi(x) + A 2 y2(x) given 
in (T23]) or ff24|) is well defined function on the interval [0, d] where d = CR 2 and R is 
the stellar radius, the quantities u, A, p, p and E are nonsingular and continuous. If 
K > l(a > 0) or K < — l(a < 0), then it is clear from ( \22ch that p remains positive 
in the region 

< 2{K 2 - 1) 



(3 + x) 

for positive constant C, which restricts the size of the configuration. We require that 
the pressure must vanish across the boundary r = R which implies that 



(K + 1 + CR 2 ) 
: {K+iy(l + CR 2 ) 



1 - K 2 + CR 2 aKCR 2 

__ CR2 + (K + iy(i + cR 2 ) + 2(K + iy(i + cR 2 y ~ 



where y is given by fT23l or ff24^) . Essentially this places a restriction on the constants 
Ai and A 2 . The interior metric ([1]) must match to the exterior Reissner-Nordstrom 
line element 

( 2M Q 2 \,o / 2M Q 2 \ _1 l9 9/ „ 9/1 ,,9v 
ds 2 = - 1 + -V dt 2 + 1 + -V dr 2 + r 2 (dd 2 + sin 2 ^# 2 
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at the boundary r 



R. This requirement implies that 



1 - 



2M Q 2 



A 2 [A iyi (CR 2 )+A 2 y 2 (CR 2 )} 2 



( 



2M Q 2 



) 



-i 



1 + CR 2 



1 



(1 + kCR 2 ) 2 



This gives the relationships between the constants A lt A 2 , k(or K), A and C. We 



to ensure the continuity of the electric field intensity across the boundary. This shows 
that continuity of the metric coefficients and matter variables across the boundary 
of the star is easily achieved. The matching condition at the boundary may place 
restrictions on the metric coefficients v and its first derivative for uncharged matter; 
and the pressure may be nonzero if there is a surface layer of charge. However there are 
sufficient free parameters to satisfy the necessary condition that arises from a particular 
physical model under consideration. 

We are in a position to investigate the gravitational behavior of this model in the 
interior of the star for particular choices of the parameter values in (126|) . The behaviour 
of the stellar model is illustrated in terms of graphs of the matter variables and the 
gravitational potentials. We have generated these graphs with the assistance of the 
software package Mathematica. For simplicity we make the choices K = 5, A = C = 1, 
di = a 2 — 1 and a = 1, over the interval < r < 1, to generate the relevant plots. 
In Fig. 1 and Fig. 2 we have plotted the metric functions e 2v and e 2A , respectively. 
It can easily be seen tha the gravitational potentials remain regular in the interior of 
the star for < r < 1. In Fig. 3 we have the behaviour of the energy density p, and 
Fig. 4 gives the representation for the isotropic pressure p. We observe that the energy 
density and the pressure are positive and monotonically decreasing functions in the 
interior of the star. The electric field intensity E 2 is given in Fig. 5 which is positive 
and monotonically increasing. Thus the quantities p, p, E, e 2u and e 2A are continuous, 
regular and well behaved throughout the interior of the star. In Fig. 6 we have plotted 
■j- on the interval 1 < r < 1. It can be observed from Fig. 6 that the speed of sound 
is always less than unity. Consequently the speed of the speed of sound is always less 
than the speed of light and causality is not violated. Therefore we have demonstrated 
that there exist particular values for the parameters so that the solution (1261) satisfies 
the requirements for a physically reasonable charged star. 



must have 



Q 2 (R) 



aKC 2 R 6 



(K+1) 2 (1 + CR 2 ) 2 
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8 Discussion 



We have found new solutions (fill) to the Einstein-Maxwell system (jlj), by utilising the 
coordinate transformation (171) , that do not have an uncharged analogue. These solu- 
tions are given in terms of elementary functions; other solutions are possible in terms of 
a general series. Consequently other new exact solutions (I22p to the Einstein-Maxwell 
field equations were found in terms of special functions, namely hypergeometric func- 
tions. The electromagnetic field may vanish in the general series solutions and we 
can regain uncharged solutions. It is possible for hypergeometric functions to be ex- 
pressed in terms of elementary functions for particular parameter values. We used this 
feature to find two classes of exact solutions (12"5|) and (1241) to the Einstein-Maxwell 
system in terms of polynomials and product of polynomials and algebraic functions. 
The simple form of the solutions found facilitate the analysis of the physical features 
of a charged sphere. For particular parameter values we showed that it is possible to 
model a physically acceptable charged relativistic sphere. 

We should emphasise that the solutions found in the paper depend crucially on the 
transformation ([7]) in which k ^ and k ^ 1. Consequently we cannot regain the 
Schwarzchild interior metric (k = 1) or the family of metrics of Hansraj and Maharaj 



19| (k = 0). A different coordinate transformation from (jTJ, allowing for k = and 
k — 1, must be utilised to regain previously known solutions; a paper outlining this 
further new class of Einstein-Maxwell solutions is under preparation. Clearly such 
solutions are possible as the following example illustrates. For the choice of metric 
function we can show that the system (jl]) admits the particular exact solution 

^ - OTW (28a) 

e 2v = 1 (28b) 



C[6(l - 2k) + x{l -2k- Ilk 2 ) - 6k 



' = — 2(l + s) a (28C) 

C[2(2k - 1) + x(2k + 3k 2 - 1) + 2k 2 x 2 } , ooJ . 

p = WTx) 2 (28d) 



, C(l - k) 2 x 

E ~ 'i + xY (28e) 



When k = 1 then E = and we have uncharged matter with the line element 

ds 2 = -dt 2 + 1 -^dr 2 + r 2 (d9 2 + sin 2 dd<j) 2 ) (29) 

1 + Cr 
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with equation of state p+3p = 0. Thus we have regained the familiar Einstein universe. 
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Figure 4: Isotropic pressure p 




0.2 0.4 



Figure 5: Electric field intensity E 2 
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